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Abstract 

We extend the Doering-Constantin approach to upper bounds on energy 
dissipation in turbulent flows by introducing a balance parameter into the 
variational principle. This parameter governs the relative weight of differ- 
ent contributions to the dissipation rate. Its optimization leads to improved 
bounds without entailing additional technical difficulties. For plane Couette 
flow, the high-i?e-bounds obtainable with one-dimensional background flows 
are methodically lowered by a factor of 27/32. 
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I. INTRODUCTION 



While it is not feasible to obtain exact solutions to the equations of motion for turbulent 
flows, it is possible to derive mathematically rigorous upper bounds on certain quantities 
characterizing turbulent flow fields, such as the rate of energy dissipation in shear flows or the 
rate of heat transport in turbulent convection. A detailed theory of upper bounds has been 
developed by Howard |jj and Busse The bounds obtained by their method, although 

still considerably higher than experimentally observed values, have until now resisted any 
attempt of further improvement (see, e.g., 0). 

Recently a quite different approach has sparked renewed interest in the theory of upper 
bounds. Instead of starting from the usual Reynolds decomposition, Doering and Con- 
stantin |||7j utilized an idea put forward already by Hopf |§ and decomposed a turbulent 
flow field into a stationary "background flow" and a fluctuating component. The background 
flow should be regarded as an arbitrary mathematical auxiliary field that merely has to carry 
the boundary conditions of the physical flow field, rather than as a time average of the actual 
flow. This method turns out to be fairly versatile [|?T,|9|-|T2|| . A given background flow imme- 
diately yields an upper bound, provided it satisfies a certain spectral constraint. The best 
bounds obtainable in this way can therefore be computed from a variational principle for 
the background flow 0. A formal connection between this approach and the Howard-Busse 
theory has been elucidated by Kerswell [JIB] . 

The objective of the present paper is to point out that the variational principle suggested 
by Doering and Constantin can still be improved, without introducing additional technical 
complications. Although our arguments are of a more general kind, we restrict the discussion 
to the concrete example of energy dissipation in plane Couette flow. In this case the benefits 
of our formulation of the variational principle are particularly obvious: the laminar flow now 
becomes an admissible background flow for Reynolds numbers up to the energy stability 
border, which is just what should be required on intuitive grounds. This eliminates a 
shortcoming of the original Doering-Constantin approach. For asymptotically high Reynolds 
numbers, and employing one-dimensional background flows only, the improved principle 
yields upper bounds on the rate of energy dissipation that are systematically by a factor of 
27/32 lower than those that can be computed from the original principle. 



II. FORMULATION OF THE VARIATIONAL PRINCIPLE 

We consider an incompressible fluid confined between two infinitely extended rigid plates. 
The lower plate at z = is at rest, whereas the upper one at z — h moves with constant 
velocity U in the positive x-direction. The velocity field u(x, t) satisfies the Navier-Stokes 
equations 

<9<u + u- Vu + Vp = uAu (1) 

with 

V • u = , (2) 

where p is the kinematic pressure and v the kinematic viscosity. We impose no-slip boundary 
conditions (b.c.) on u at z = and z = h, 
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u(x, y, 0, t) = , u(x, y, h,t) = U~x 



(3) 



(x is the unit vector in x-direction), and periodic b.c. on u and p in x- and y-direction; the 
periodicity lengths are L x and L y . The time- averaged energy dissipation rate per mass is 
given by 



£r = <HVu| 2 )) T = i^i 



V 

ttJa 



(4) 



We employ the notation (•) = Ind 3 x (•) for the volume average and (-) T = i J^dt (•) for the 
time average, Q = L x L y h denotes the periodicity volume, and u^j symbolizes djUi(x,y, z,t) 
with i,j = x,y,z; summation over repeated indices is implied. The goal is to formulate a 
variational principle for rigorous upper bounds on the long time limit of Et, 

e = limsup Et , (5) 

T^oo 



or on the non-dimensionalized quantity 

U 3 h 



Ce(Re) = TT^TT . ( 6 ) 



where Re = Uh/v is the Reynolds number. 

The approach by Doering and Constantin ||[7| rests on a decomposition of the velocity 
field u(x, t) into a background flow and fluctuations, 

u(x,t) = U(x) + v(x,t) . (7) 

The background flow U is a stationary and divergence-free vector field satisfying the physical 
b.c: U(x, y, 0) = 0, U(x, y, h) = Ux, and U is assumed to be periodic in x- and ^-direction; 
otherwise it is completely arbitrary. For the divergence-free fluctuations v we then have 
homogeneous b.c. for all instants t > 0, i.e., v(x, y, 0, t) = v(x, y, h, t) = 0, and v is periodic 
in x- and ^/-direction. With the usual manipulations one arrives at the following equations 
for u and v: 

|Vu| 2 ) = z/(|VU| 2 ) + 2z/(%^ b -) + z/(|Vv| 2 ) , (8) 



d 
di 



v 2 ) + (v(U- VU)) + (v ■ (VU) 



/sym 



— V 



Vvl 



(9) 



where 



(VU) 



sym 



\Uj\i + Ui\jj/2. At this point, Doering and Constantin use eq. (H) 

to eliminate the cross background-fluctuation term in eq. @, namely, they form [eq. ([§])] + 
2 • [eq. (S)]. It is crucial to note that this implies putting a certain fixed weight on the 
different contributions to the resulting expression for Et- In order to avoid such a weighting 
we introduce a new degree of freedom and consider 



[eq. (8)] + a ■ [eq. 



a > 1 



(10) 
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with the balance parameter a. Apart from this modification, we now follow closely the spirit 
of the original Doering-Constantin approach, which can always be recovered by setting a = 2. 
Equation (|T0|) leads to the energy balance equation 



£T= ,(|VU| 2 )-|,(iv(,T) 2 ) + | i (iv(,0f 



a — I 



HVv| 2 + v(VU) sym -v + f -v 



with 



a - 2 

f = U • VU i/AU 



(11) 



(12) 



which becomes the inequality (see (||)) 
e < z/(|VU| 2> 
— a lim 



lim inf / — f 

T^oo \ Q, Jq 



d 6 x 



a-l 



z/|Vv| 2 + v(VU) sym -v + f -v 



T 



(13) 



For the evaluation of (|T3|) we have to distinguish two different cases. 
Case a) f is a gradient 



In this case the term linear in v on the rhs. of ([13]) vanishes and the resulting bound on the 
energy dissipation rate is given by 



e< inf |z/(|VU| 2 \} 
" u,«>i I V / i 



(14) 



provided U complies with the following constraints: 

i) U is a divergence-free vector field satisfying the physical b.c. and f is a gradient, 

ii) the functional 



#u,a{w} = - / d 3 x 



a v |Vw| 2 + w ■ (VU) sym • w 



(15) 



is positive semi-definite, -£/u,a{w} > 0, for all stationary divergence-free vector fields w 
satisfying the homogeneous b.c. 

The condition ii is equivalent to the statement that all eigenvalues A of the hermitian eigen- 
value problem 



a - 1 



A V = -2 — z/AV + 2 (VU) sym • V + VP 



= V ■ V , V satisfies the homogeneous b.c. , 



(16) 



are non-negative; P is a Lagrange multiplier for the divergence condition. Following Doering 
and Constantin 0, who consider the eigenvalue problem without an adjustable parameter 
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(a = 2), we denote the requirement that all eigenvalues of the eigenvalue problem ([TJ) be 
non-negative as 

spectral constraint: all A > . (17) 



Case b) f is no gradient 

In this case we bound the second term on the rhs. of ([T^) by 
'a — 1 



i I d 3 x 



Q Jn 



HVvr + v(VU) sym -v + f -v 

a 



> inf { - I d 6 x 



1 2 

- v | Vw| + w • (VU) sym • w + f • w 



(18) 



where we seek the infimum in the space of all stationary, divergence-free vector fields w 
satisfying the homogeneous b.c. If V is the minimizing vector field, it solves the Euler- 
Lagrange equations 







z/AV + 2 (VU) sym • V + VP + f , 



= V • V , V satisfies the homogeneous b.c. 
so that the inequality (|T8"D becomes 



^-I^|Vv| 2 +v(VU) sym -v + f -v 



> ^- [ d 3 x f ■ V 

T ~ 2 Jn 



The requirement 

iJu.o{w} > for all divergence-free w/0 satisfying the homogeneous b.c. 



(19) 



(20) 



(21) 



guarantees the uniqueness of the solution V to the Euler-Lagrange equations ( |T9"D and en- 
sures that this solution indeed minimizes the rhs. of ([18]). The uniqueness can be shown as 
follows: if there were any other solution W ^ V, then the difference W — V would be a 
nonvanishing eigensolution to the eigenvalue problem (|1^) with eigenvalue A = 0. This con- 
tradicts the requirement (|2~T|), which implies that all eigenvalues of ( |TED are strictly positive. 
The minimizing character of V stems from the convexity of the expression in curly brackets 
on the rhs. of (|18D . 

Putting all things together, we see that in this case b) the resulting bound on the energy 
dissipation rate is given by 



e< inf lu(\VU\ 

U, a>l 



^(V.(U.VU)) + ^,(V.AU)} 



(22) 



provided 

i) U is a divergence-free vector field satisfying the physical b.c. and f is no gradient, 

ii) V is a solution to the Euler-Lagrange equations ([19|) , 
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iii) the functional CP~5|) is strictly positive definite, Hu ta {w} > for all divergence-free 
vector fields w ^ satisfying the homogeneous b.c. Equivalently, all eigenvalues A of 
the eigenvalue problem (|T6|) must be positive, so that now we have the 

spectral constraint: all A > . (23) 



The main technical problem encountered in the original Doering-Constantin approach is to 
verify their spectral constraint for a given background flow, i.e., to show that all eigenvalues 
of the eigenvalue problem flTE| ) with a = 2 are non- negative (case a) or positive (case b). In 
this respect our formulation generates no additional complications because our constraining 
eigenvalue problem is formally identical with the original one, the only difference being that 
the kinematic viscosity is rescaled by a factor 2 (a — I) /a. Since —A is a positive definite, 
unbounded operator, the spectral constraint enforces the positivity of this scaling factor, as 
anticipated in eq. ([U]) by the condition a > 1. 

It should be noted that our formulation differs from the original one in two major points. 
Firstly, if AU is no gradient then the minimizing field V must be determined by solving the 
Euler-Lagrange equations (|l^) even if U • VU is a gradient, since the contribution to f that 
is proportional to AU vanishes only if a = 2, see eq. (|T2"|). Secondly, our variational principle 
concerns not only the background flow U but also the balance parameter a. We will see 
that in the case of one-dimensional background flows both issues can be dealt with easily; 
solving the Euler-Lagrange equations and minimizing over a will be independent from the 
intricate problem of verifying the spectral constraint. Hence, the additional freedom gained 
by the balance parameter a entails no additional difficulties, but will result in an improved 
bound on e. 



III. PLANE COUETTE FLOW WITH ID BACKGROUND FLOWS 

To become more specific we restrict the following discussion to one-dimensional back- 
ground flows, i.e., to flows which can be described by a profile function 0, 

U = Ety(C)i; 0(0) = 0, 0(1) = 1, (24) 

depending on the dimensionless variable ( = z/h. Clearly, such a U is a divergence-free 
vector field satisfying the physical b.c. Additionally we impose 0(C) = 1 — 0(1 — C)> so that 
the profile is adapted to the symmetry of the physical problem. 

a) Laminar profile, 0(£) = C, 

In this case U is the laminar solution to the equations of motion (|]) and (0), and we have 
f = 0. Hence, case a) of our variational principle applies. Provided all eigenvalues of the 
eigenvalue problem ([16]) are non-negative, the estimate ([14]) yields 

U 2 

e = v — or c £ = Re 1 , (25) 
h 

respectively. We have strict equalities here since upper and lower bounds on the long time 
limit of Et coincide, as long as the laminar profile fulfills the spectral constraint J7[. We are 
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thus led to the following question: what is the maximal Reynolds number Re (or inverse 
kinematic viscosity up to which the laminar profile is admitted as a valid testprofile for 
the variational principle? 

Since the viscosity enters into the eigenvalue problem ( |l6l) only in rescaled form, we 
can tune the Reynolds number by suitably adjusting the value of a. In order to find the 
maximal "critical" Reynolds number for the laminar profile we have to set a = oo, so that 
the scaling factor 2 (a — l)/a of the viscosity takes on the highest value possible, namely 2. 
The eigenvalue problem then becomes 



AV = -2z/AV + - 
h 



/0 1 
I ■ V + VP 



V i o o 

= V ■ V , V satisfies the homogeneous b.c. (26) 

This is exactly the eigenvalue problem appearing in energy stability theory for the plane Cou- 
ette flow |TJ],rj|. Hence, the laminar profile is an admissible testprofile up to the Reynolds 



number PeEs characterizing the energy stability border fl4| . |15l 



2Vl707.77 ~ 82.65, (27) 



and consequently (^) is valid for all Re < Pess- 

In the original Doering-Constantin approach (a = 2) one had to discard the laminar 
profile already for Reynolds numbers Re > PeEs/2, which necessarily led to non-optimal 
upper bounds on c e for Reynolds numbers PeEs/2 < Pe < PeEs- The introduction of the 
balance parameter cures this obvious shortcoming and guarantees that c £ as obtained from 
the variational principle has the known 1/ .Re-behaviour up to Re = Re^s- 

b) Non-laminar profile, 0(£) ^ £ 

In this case we have 

f = -— i/|U"(C)* with f(()^0. (28) 
a h l 



The boundary conditions (|24|) together with the symmetry condition ensure that f cannot 
be a gradient if a ^ 2, so that now we have to resort to case b) of our variational principle. 
Because of the special form that the inhomogeneous term f acquires for non-laminar profiles, 
one can find an analytical solution to the Euler-Lagrange equations fliD|): 



v = ^^£/[C-<KC)]x, (29) 



2 a- 1 



l Jo 



(30) 



Note that V vanishes if a = 2, so that we do not need to distinguish between the cases a^2 
and a = 2. Moreover, the rhs. of ( p9|) is proportional to the difference between <f)(() and 
the the laminar profile (, which leads us back to the laminar case a) in the limit — > (. 
Provided all eigenvalues of ([TB]) are positive, the estimate (|22|) yields 
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£ < inf 

<j>, a>l 



4(a-l 



D{0} 



the inequality for c e reads 



c £ < inf 

<l>, a>l 



1 + 



4fa-l 



D{0} 



LP 
h 2 



Re' 1 



(31) 



(32) 



Here we have employed the abbreviation D{<p} for the functional 



D{0} = [\ [0'(C)f - 1 
Jo 



(33) 



This functional is strictly positive for <p(() ^ (, which means that for each non-laminar 
profile and each a > 1 the factor 1 + ^-i) ^i^} appearing in ([Jl]) and fl32|) exceeds one. 
Therefore, each non-laminar profile produces a bound that is strictly higher than the laminar 
bound (0). 

In a manner analogous to the procedure for the laminar profile one has to investigate 
on the basis of the eigenvalue problem (|16|) up to which Reynolds number a given profile 
is admissible as a testprofile. Because of the positive definiteness of —A, the spectrum of 
the hermitian operator defined by the rhs. of fli~6|) is lowered when the rescaled kinematic 
viscosity is decreased. Thus, for a given we define the critical Reynolds number R c {(f)} as 
that Reynolds number where the lowest eigenvalue A of ([16|) with fixed balance parameter 
a = oo (i.e., 2 (a — l)/a = 2 is maximal) passes through zero. For Re > R c {4>} one has 
to discard as a testprofile for the variation. In addition, if a > 1 is finite, one finds the 
constraint 



Re < - — - R c {(p} 



(34) 



The factor a 2 /(a — 1) in ([31]) and (|32]) has a local minimum for a = 2 and increases 
monotonically with a for a > 2. A given thus produces the following upper bound on c e : 



c, < 



1 + 



4(a r 



D{0} 



Re' 1 for < Re < i? c {0} 



with 



hence 



c e < 



R c {<l>}-Re 



1 + D{4>}} Re' 



for < Re < \R c {(p} 

for \R c {<p} <Re< i? c {0} 



1 + 



D{4>}RM} 



4(R c {(j>}-Re)Re 



Re~ 



for < Re < ±i? c {0} 

for \R c {<p} <Re< R c {(p} 



(35) 



(36) 



(37) 



In this way we have accomplished the optimization of the balance parameter and are left 
with the task of varying the profile function. But it is possible to deduce some general 
statements even without solving the variational principle for 0. To this end, we denote the 
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expression on the rhs. of ([37]) as c e (Re). This is a continuous function of the Reynolds 
number, and even continuously different iable at Re = R c {4>}/2. For every given profile, 
c £ (Re) has exactly one local minimum in the whole interval < Re < R c {<p}; this minimum 
appears in the upper half interval R c {4>}/2 < Re < R c {<p}. Because the variational principle 
tests all profile functions <fi satisfying the required conditions, the minimum point determined 
by a particular is the only point that this profile could possibly contribute to the resulting 
upper bound on c £ (Re). Thus, each leads to a point in the (Re, c e )-plane, 

(f> — > (iW0},c e (iW0})) • (38) 
The Reynolds number R m i n {4>} of the minimum point can be expressed as 

i? min {0} = x o {0} J R c {0} , (39) 
where £o{0} is the unique (real) zero Xq with 1/2 < x < 1 of the cubic polynomial 

x 3 -2x 2 + (l- £>{</>}) x + ^£>{0} = . (40) 

This follows directly by minimizing c £ (Re) with respect to Re. Although an analytical 
expression for the desired zero of (|40"D is available for all D{<p} > 0, we will discuss only two 
limiting cases. In the laminar limit <p(C) — > £ the functional D{<p} vanishes and the zero 
xo{(p} tends to 1. Inserting i? m i n {0} into ( |36|) shows that a m \ n tends to infinity in this limit, 

lim a min = oo . (41) 

Re \ Re es 

Therefore, the optimal bounds on c £ resulting from the variational principle will be con- 
tinuous at the energy stability border. On the other hand, for sufficiently high Reynolds 
numbers the spectral constraint ( [23] ) singles out those profile functions as admissible for the 
variational principle that have large slopes within thin boundary layers and remain almost 
constant in the interior [7J. Hence the functional -D{0} tends to infinity in the asymptotic 
limit Re — * oo, which implies that the zero xo{4>} tends to 2/3 and a m i n approaches its 
asymptotic value 3, 

lim a min = 3 . (42) 

The optimal value of the balance parameter a, considered as a function of the Reynolds 
number, is monotonically decreasing from infinity at the energy stability border -ReES to its 
asymptotic limit a m = 3 and never reaches the value a = 2 that has tacitly been taken for 
granted before. 

In the original approach a given profile <ft yields an upper bound on c £ in a similar 
way, but the range of Reynolds numbers within which this profile is admissible is only half 
as large as in our case. Since, as remarked before, the factor a 2 /(a — 1) appearing in (|32|) 
has a minimum for a = 2, in the common interval this bound agrees with ours: 

c £ < c £ (Re) for < Re < - R c {0} . (43) 
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Here the minimizing point is (R c {<p}/2, c £ (R c {4>}/2)), which has to be compared with (j38[). 
Therefore, in our formulation a given profile function leads to a lowered bound on c £ at 
a higher Reynolds number Re. The improvement due to the balance parameter a relies on 
the fact that each <fi becomes admissible within an enlarged .Re-interval. Note that in both 
cases exactly the same eigenvalue problem must be solved to check the spectral constraint. 

Under the assumption that in both formulations the upper bounds on c e are approaching 
constant values in the asymptotic limit Re — > oo, dubbed c £)OC and c e D ^, respectively, we 
can calculate the relative factor 

„ - | itn CeM^jRcW) _ 27 

^-c-dc-^ c - £(jRc{0} /2) -32" (44j 

Considering the same class of one-dimensional background flows, our formulation yields 
asymptotic upper bounds on the rate of energy dissipation that are systematically lowered 
by a factor of 27/32 compared to those that can be computed from the original Doering- 
Constantin approach. 



IV. ILLUSTRATION 



In this section we wish to illustrate the general statements derived in the previous section. 
Our aim is not to calculate the best bounds that the improved variational principle has to 
offer, but rather to elucidate the new aspects in a simple way. Instead of verifying the 
spectral constraint for a given on the basis of the eigenvalue problem (|T6|), one can check 
the positivity of the functional ifu,a{w}, 



#u,a{w} = - / d 3 x 
\l Jn 



1 IV7 \*m U J 
— V VW +—(f>W x W z 

h 



(45) 



by means of the inequality 

— / a x — <p w x w z 
il Jn h 



< 



Vwl 



h/2 {J 

dzz — 
h 



(46) 



This estimate can be shown by repeatedly using Schwarz' inequality and utilizing the sym- 
metry property of the profile function |Tl| (see also refs. ||0]). Thus, 



#U,a{w} > 



Uh 



V 



,1/2 




/ dCC ^ 


(|VW| 
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(47) 



If the profile function <f>(C) satisfies the condition 



1/2 



dCC 1^(01 < 2 ^2- 



Re 



-1 



or 



[ 1/2 d((\<P'(0\<2V2 — Re- 1 
Jo a 



(48) 



then the non-negativity (or the positivity) of the functional i?u,a{w} is guaranteed. It 
should be realized that (|4"8"D is a sufficient but not a necessary condition; it is more restrictive 
than the spectral constraint. In addition to the b.c. (|24]) and the symmetry requirement we 
assume </>'(C) ^ 0> so that the modulus signs || in (TO) can be skipped. 
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After the replacement of the variational principle's spectral constraint (condition ii in 
case a) and condition iii in case b), respectively) by the sharpened profile constraint (fj5|), a 
certain can, in general, no longer be admitted as a testprofile for Reynolds numbers up to 
R c {4>} as defined in the previous section. Rather, ( [Hf) leads to the border 

Jo rf CC0 (0 

For example, the laminar profile <f)(Q = £ yields 

Rl{(} = 16^2 ~ 22.63, (50) 



which has to be contrasted to (p7|), i.e., -R C {C} = R e ES — 82.65. Thus -Rf{C} is by about 
a factor 4 smaller than R c {(}. Nevertheless, it is instructive to consider the variational 
principle with the stronger constraint fl4~8|) since then even the variation of the profile <fr can 
be done analytically. For finite a > 1 we now have the constraint 

Re < R s c {(j)} (51) 

instead of (|34]), which is exactly ([48]). 

Without going into the cumbersome technical details (see also ref. [T(J for comparison) 
we summarize our results. We have to distinguish three .Re-ranges, 

(I) < Re < 16V2 , 

(II) 16V2 < Re < 20V2 , (52) 

(III) 20V2 < Re < 00 . 

In these ranges the minimizing profile function behaves as sketched in Fig. 1: in (I) <p is 
equal to the laminar profile, <p{Q = (. In (II) the profile develops parabolic boundary layers 
which reach the middle of the domain Q (i.e., £ = 1/2), while the slope at £ = 1/2 decreases 
from one to zero with increasing Re. When Re is increased beyond 20v2, the thickness S 
of these parabolic boundary layers is getting smaller and smaller; in the interior remains 
constant, = 1/2 for 5 < ( < 1 — 5. Asymptotically, 5 vanishes as ~ 1/ Re. The analytic 
expression for 5 in (III) reads 



5 = 4y/2 3^JRe -2V2 + ^Re - 18y^ 8^/2 
Re jRe-2V2 + jRe-18V2 Re 



Correspondingly, the functional D{<p} defined by eq. ( |33"D diverges for Re —> 00, as antici- 
pated in section 3. 

The minimizing balance parameter a turns out to be 



00 



(i-e-, 2 *=! = 2) (I) 



16V2 

Re-lQy/2 



(H) . (54) 
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For Re > 16^/2 this parameter decreases monotonically from oo to its asymptotic value 
depicted in Fig. 2. 
Finally, the resulting bounds on c £ are given by 



l<c e (Re)< 



35 « 

i - £ (n) 



3 



Re- 2V2+^J Re-l&V2 



+ 7 . 1 . =v» (I") 



(55) 



16 V^ 3 Vfle-2>/2+Vfle-18v^ (V^- 2^2+ V-Re- 18^2 



Note that the upper bound is continuous but not continuously differentiable at Re = 16\/2, 
see Fig. 3. The limit Re — > oo yields 

^°°5< J & ^(ite) = ^- (56) 

As before, we denote the upper bound on c e by c e (Re). 

These results have to be compared with those derived in ref. [|Kj], where the sharpened 



profile constraint ( (48|) was used to calculate bounds on c £ analytically within the framework 
established by Doering and Constantin (a = 2 fixed). The analogous three .Re-ranges found 
there are 

(I') < Re < 8V2 , 

(II') 8^ < Re < 12V2 , (57) 
(III') 12^2 < Re < 00 ; 

in particular, the first range (I') is only half as wide as in our case. The corresponding 
bounds on c £ are given by 

^ < c £ (Re) <{ + (HO • (58) 



(III') 



9V2 



The improvement due to the parameter a is measured by the ratio of the rhs. of (pq ) and 
the rhs. of (1581), 



^ s WW) • < 59) 

This ratio has a minimum at Re = 16y/2, 

/ r-\ 9\/2 9 

§ ( 16VI ) = wv=2 = is * a56 • (60) 

and increases monotonically with increasing Re to its asymptotic value 27/32, 

lim g{Re) = = — ~ 0.84 . (61) 

Re^oo yK ' 32V2 32 v ; 

A graphical comparison of the bounds fl5"5]) and (p8|) is shown in Fig. 3. 
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V. CONCLUSION: BACKGROUND FLOW AND BALANCE PARAMETER 

BELONG TOGETHER 



In this paper we have been concerned with an improved formulation of the Doering- 
Constantin variational principle for bounds on turbulent energy dissipation, rather than 
with its solution. The necessity of devoting the utmost care to an optimal formulation is 
obvious: what has been lost by a non-optimal formulation of a variational principle can not 
be regained by even the most sophisticated techniques for solving it. 

The feasibility of improving the background flow method rests on two key observations. 
The first of these is that if one refrains from what appears to be straightforward, namely, 
if one does not eliminate the unwanted cross background-fluctuation term from eq. @ but 
rather keeps this term with a certain weight quantified by the balance parameter a, one 
gains a new freedom that can be exploited to improve the bounds. The second observation 
is that the required solution to the Euler-Lagrange equations (|B^) can easily be written 
down analytically if one restricts oneself to background flows which only have a height- 
dependent profile but no spanwise structure. Then the technical difficulties encountered in 
our formulation are the same as those of the original approach — in both cases one has to 
solve the same eigenvalue problem in order to check whether a given profile is an admissible 
candidate for the computation of the bound, albeit in our case the kinematic viscosity is 
rescaled — , but the improved formulation leads to systematically lowered bounds. 

Finally, we wish to stress that the usefulness of introducing the balance parameter is 
by no means restricted to the case of plane Couette flow. It should be obvious by now 
that the arguments employed in sections 2 and 3 can easily be adapted to other problems. 
When computing upper bounds on any property of a turbulent flow that is amenable to 
the background flow method, background flow and balance parameter should henceforth be 
regarded as belonging together. 
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FIGURES 



FIG. 1. Background flow profiles resulting from the variational principle with the sharpened 
profile constraint (f48|). Solid line: < Re < 16\/2; short dashes: Re = 20^2, long dashes: 
Re = 64^/2. 



FIG. 2. Balance parameter a resulting from the variational principle with the sharpened profile 
constraint ((4S|). 



FIG. 3. Bounds on the dimensionless energy dissipation rate c £ (Re) derived from the variational 
principle with the sharpened profile constraint (|48|). Solid line: lower bound, 1/Re; long dashes: 
upper bound with optimized parameter a (see inequality (|55|)); short dashes: upper bound obtained 
in ref. [ 10 1 for a = 2 fixed (see inequality 
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